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I. INTRODUCTION 
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Transport coefficients measure the efficiency with which conserved quantities are transported through a medium, 
over distances that are long compared to the microscopic relaxation scales of the system (for a review see [1]). Direct 
applications include the early universe and the quark gluon plasma. The calculation of transport coefficients is 
also important from a purely theoretical standpoint. They characterize linear deviations from equilibrium, but are 
calculated using the familiar methods of equilibrium field theory. Results could therefore be used as a check on 
Qh| calculations based on purely non-equilibrium methods. 

The calculation of transport coefficients in gauge theories is hugely complicated by the occurrance of pinch- and 
collinear-singularities [3,[1,11|. These singularities can be regulated by using hard thermal loop propagators. However, 
this remedy produces infinite sets of graphs which all contribute at the same order. In order to include all leading 
order contributions, these infinite sets of graphs need to be resummed. The problem is to develop a technique to 
perform the resummation that avoids double counting and respects all of the symmetries of the original theory. 

The complete leading order calculation of electrical conductivity and shear viscosity, in both QED and QCD, was 
done in [||. This calculation is not obtained directly from quantum field theory, but is derived from kinetic theory. It 
is of interest to understand the connection between the kinetic theory approach, and a calculation based on quantum 
field theory. One motivation is that quantum field theory might provide a better framework than kinetic theory for 
calculations beyond leading order. The equivalence of the two approaches has been demonstrated for scalar theories 
[H, & 0, HI • In QED, field theory based calculations have been done usin g a d irect ladder summation (UHijj dynamical 
remormalization group methods and other diagrammatic methods [12L fl3| . A large Nf leading log calculation of 
conductivity and shear viscosity was done in 11411 using the 2PI effective theory. The conductivity has been obtained 
at leading order from the 3PI effective action [IE fl6| . In QCD, the field theoretic calculation of shear viscosity has 
only been done at leading log order [H, [It}. In this paper we present the first calculation of the full leading order 
QCD shear viscosity using quantum field theory methods. We show that the calculation can be organized naturally 
using the 3PI effective action. The result provides strong support for the use of nPI effective theories as a method to 
study the equilibration of quantum fields. 

The n-PI effective theory, in which the n-point functions are treated as variational parameters, is a natural method 
to organize the calculation of transport coefficients. In general, a consistent resummation requires the computation 
of the nPI effective action for infinite n. However, there is an equivalence heirarchy that simplifies the structure of 
the calculation [18]: to 2-loop order, the infinite-PI effective action is equivalent to the 2PI effective action, to 3-loop 
order, the infinite-PI effective action is equivalent to the 3PI effective action, etc. 

Truncations of the effective action produce problems with gauge invariance. Even though the effective action is 
consistent with the global symmetries of the theory, the Ward identities associated with the gauge symmetry may 
not be satisfied for the self consistently determined vertex functions [HI, Hoj. To address this problem, we use the 
resummed effective action, which is defined with respect to the self-consistent solutions of the n-point functions 
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There are three different types of n-point functions involved in the calculation: (1) the self-consistent solutions of 
the equations of motion; (2) 'mixed' n-point functions obtained by taking functional derivatives of the self-consistent 
solutions with respect to field expectation values and; (3) 'external' n-point functions obtained by taking functional 
derivatives of the resummed effective action with respect to field expectation values. For QED, it has been shown that 
the 'external' n-point functions satisfy the usual Ward identities [H, IHj]. At the exact level, all of these definitions 
are equivalent to each other. Integral equations for type (1) vertices are obtained from the equations of motion of 
the effective action (see Eqn. (fT7|) ). Integral equations for type (2) vertex functions are obtained by functionally 
differentiating the equations of motion with respect to field expectation values (see, for example, Eqn. (|32|) ). 

This paper is organized as follows. In section [Til we define some notation. In section [Til] we derive an expression for 
the shear viscosity in terms of the integrand that gives the gluon polarization tensor. In section [IV] we give the 3PI 
effective action to 3-loop order. In section [V] we define some 'external' n-point functions and obtain an expression for 
the 'external' 2-point function. In sections PvTl and IVDl we derive integral equations for the relevant 'mixed' vertices, 
and the self consistent vertices. In section IVIIII we show that the kernels of these equations can be written as the 
square of the sum of the amplitudes that correspond to all physical scattering and production processes. In section 
HXl we present our conclusions. In Appendix A we define the notation used in the Keldysh representation of finite 
temperature field theory. In Appendix B we show that the expansion of the 'external' polarization tensor produced 
by the 3PI formalism contains all of the terms that would be produced by a Wick expansion. In Appendix C we give 
some details of the calculation presented in section IVIIII 



II. NOTATION 



We use: 



# = diag(l, -1,-1,-1), (1) 

N B (po) = 1 + 2n b (p ) , N F (p ) = 1 - 2n f (p ) , 
/ dP:= /(2^' / := /2^(L)3' 



v 3 y p z 

We work in the Feynman gauge and write the QCD Lagrangian as: 

C = -\F« V F^ a - \{d*Al? + i^D^ - f} a d^D^r , (2) 

F% = d^Al - d v Al +gf abc AlAl , 
D^ = {d„-igAp a )^, 

The classical action is: 

Sct[il>,$,A,ri,fj] = Jd 4 x£. (3) 

The group factor notation for SU(N) is: 

fundamental representation : C F = (TV 2 - l)/(2N) , T F = 1/2 , d F = N , (4) 
adjoint representation : Ca = N , Ta — N , d,A = N 2 — 1 . 
For simplicity we will set the coupling constant g to one throughout. 



III. SHEAR VISCOSITY 



The Kubo formula for shear viscosity is: 

j] 



1 f d OT 
20 V dqo 



r(go,0) 



qo^O 



(5) 
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P™(<Zo,0)= / dt / d 3 xe^^(t)(7r^'(t,x),7r^'(0)) 



where 7r ZJ is the traceless part of the energy momentum tensor. 

We begin by considering the contribution to the shear viscosity that corresponds to using bare propagators and 
vertices. We will call this quantity f\. We use the Keldysh representation of the real time formulation of finite 
temperature field theory. The basic method is described in Appendix A. The derivation of the corresponding expression 
for the case of the QED conductivity is given in detail in pj3|. We obtain 1 



n=-^P J dP n/(po)(l - ra,(po))Ti:((Ao)* J c ,(P + Q,Q,P)S°(P + Q) ret (A )%(P + Q,Q,P)S°(P) adv ) (6) 
+n 6 (po)(l + n b ( Po )) Q (O )^ T (-P - Q,Q, P)D° TT ,(P + Q) ret (O )^' A '(-P - Q,Q,P) D° xx ,(P) ad 

~ (®oZ(P + Q, Q, P)G°(P + Q)ret (Oo)li(P + Q, Q, P)G°(P) adv S 



9 = 
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The vertices are defined as: 



(2irY5\K + P + Q){£l y£ T {K,Q,P) = j <?x j d A y J d^e^e^e^* {A x a {z)v% uou { X )Al{y)) , (7) 
(2n) i 5\K-P-Q)(eo)Z(K,Q,P)= [ d^x [ d% [ d l z e ikz { na {z)n% ost {x)n b { V )) , 



(2n)W(K - P - Q)(A o y c j c ,(K,Q,P) = d 4 x d*y / d*ze-**e-™e*'(M*)*^(x)Mv)) , 



ghost V 
ij 

quark V 

where {7Tq Uark , ^ghosf ^giuoni indicates the part of the traceless energy momentum tensor that is quadratic in quark, 
ghost, and gluon fields. For the ghost and quark vertices the order of the momenta is: outgoing momentum of the 
outgoing ghost/quark, incoming momentum of the gluon, incoming momentum of the incoming ghost/quark. For the 
gluon vertex all momenta are incoming. 

The integral in ([6]) contains three bubble-type diagrams: a gluon bubble, a ghost bubble, and a quark bubble. These 
diagrams contain pinching and collinear singularities that need to be resummed in order to obtain the full leading 
order contribution. In order to simplify the explanation of this point, we consider a generic bubble diagram, and draw 
the propagators as solid lines, as shown in Fig. [TJ Throughout this paper we will use stars to indicate the legs of a 
type (2) or type (3) n-point function that correspond to functional derivatives with respect to field expectation values 
(see also Fig. [6]). 




FIG. 1: A generic contribution to the shear viscosity. 



Pinch singularities are produced by the low frequency limit in the Kubo formula (Eqn. ©). When integrating a 
term of the form J dpo G ret (P)G adv (P), the integration contour is 'pinched' between poles on each side of the real 
axis, and the integral contains a divergence called a 'pinch singularity.' A set of graphs containing pinch singularities 
has the general form shown in Fig. O 




FIG. 2: Contributions to the shear viscosity from pinching singularities. 



[1] Our notation throughout this paper differs slightly from that used in [l5L[l3|. In this paper, the symbols used to represent propagators 
and vertices correspond directly to lines and dots in diagrams, with no additional factors of dzz. 
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In gauge theories, one also has collinear singularities. Fig. [3] shows a set of graphs with collinear singularities. 




FIG. 3: Contributions to the shear viscosity from collinear singularities. 

In general, pinch and collinear singularities are regulated by using hard thermal loop (HTL) resummed propagators. 
This procedure introduces extra factors of the coupling in the denominators which change the power counting. As a 
consequence, the infinite series of graphs depicted in Figs. [2] and [3] are all of the same order and need to be resummed. 
The resummation is done by solving a set of coupled integral equations that have the general form shown in Fig. [H 
When the leg on the right hand side has a star, the integral equation resums pinch singularities. The same equation 
without the star resums collinear singularities. The basic goal of this paper is to show that these integral equations 
are produced naturally by the 3PI formalism. We note that since the ghost HTL self energy is zero, the ghost diagram 
cannot be regulated in this way. However, ghosts are not physical particles and are only needed to cancel unphysical 
gluon polarizations. We will show explicitly how this works in section [Villi 




FIG. 4: An integral equation that resums singularities. 



In Feynman gauge the HTL quark, ghost and gluon propagators can be written: 

"^(P)~, -zG ab (P) = 5 -0, (8) 

Km = s ab D^(p) , - iD^(p) = _p^_i_ _ pr p^ 1 ^ - pp+s > 

P*" = - U»U V + (P" - poU^P" - poU»)/p 2 ; U = (1, 0, 0, 0) , 

P L V = ~ P T + 9^ ~ P^P V /P 2 , P£" = P^P V /P 2 . 

Dominant contributions to the shear viscosity come from hard excitations on the pinching lines, and the residue of 
the longtitudinal part of the gluon propagator is exponentially suppressed at high temperatures. Therefore, we need 
only the transverse part of the pinching gluon propagator which can be written: 

-iD^P) = ^-^-)^-. (9) 

Now we consider the vertices defined in (0. We factor the colour structure by defining: 

(fto)^ = W AT , (e )^ = ^, (Ao)« =6 cc ,Aij. (10) 
In addition, we write the bare vertices that are obtained directly from the Lagrangian as: 

(n )^b(-PAP) = f acb ^ T (-p,o,p) , (n) 
(e y acb (p,o,P) = f acb oZ(P,o,P) , 

(A ) cc ,(P,0,P) = S cc ,A (P,0,P). 

Note that we have used the same letters in and ([TT|l . In order to simplify the notation, we do not introduce 
any additional primes or tildes to distinguish the two types of vertices. The indices associated with each vertex are 
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sufficient to indicate which type of vertex is meant. Using this notation, it is straightforward to show that: 

n XT ' ij (-p,o,p) -» 2g XT P 2 P j -» p s n XsT (-p,o,p)i ij , (12) 
eg* (p, o, p) = - P 2 F = P s eg(p, o, p)/« , 

A^(p,o,p) = ^ i vy-7v, 

where the arrow indicates that the relation holds only when multiplied by the transverse projectors P^ A P^ T ' on both 
sides. We have taken the limit Q — > 0, since this produces no difficulties for factors in the numerator of the integrand. 

Using these results we can rewrite the integrand in ([6]) in terms of the integrand for the gluon self energy. We 
separate contributions to the gluon self energy from gluon, ghost and quark bubbles by writing: 

(13) 



n [C& (Q) = J dPUint[{\ s a s b (P, Q) ; ie {gluon, ghost, quark} . 
It is straightforward to show that (0 can be written: 

dPp s p s ' 



V —I o ab 
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(14) 



Nr. 



J- (nm* [gluon]*!' (P, Q) + ffint [ghost] S' (P, Q)) + ^ffint [quark]-, (P, Q) 



<f = 
go 



IV. THE 3PI FORMALISM 



At this point, we introduce a compactified notation. We use a single numerical subscript to represent all continuous 
and discrete indices. For example: a gluon field is written A\ := A?(x)\ the quark propagator is written Syi '■= 

S a p(xi,X2)', the bare 3-gluon vertex is written ^32 := (^o)aST x 3> x<2 ) e ^ c - We a ^ so use an Einstein convention in 
which a repeated index implies a sum over discrete variables and an integration over space-time variables. The free 
propagators and vertices are defined as: 

^ = • ^ = -Sk • ^ = -*£t ■ (is) 

■ 5*S cl = S(S° 12 yi o _ 5*S cl _ 5(D° 12 )-i 

132 tofoSAs&h 5A 3 ' 132 5A 2 dA 3 SA 1 5A 3 ' 

_■ PS* _ SjG^)- 1 o 6*S cl _Sn° 132 _ 5 2 (D^ 2 )~ 1 

132 SmSAsSf]! SA 3 ' 1234 6A 4 5A 3 6A 2 6A! SA 4 SA 4 5A 3 ' 

The 3PI effective action can be written (l8Ll2fij: 

T[i,,i,,A,r 1 ,fj,S,D,G,V,U,Y]=S cl [i,,i,,A,r 1 ,fj] (16) 

+ ^TvLnD^ 1 + [(r>? 2 )" 1 (D 21 - D° 21 )] - iTrLn^ 1 - iTr [(S^)" 1 (S 21 - 

- iTr LnGr 2 X - [(G? 2 ) _1 (G2i - G° 21 )] + T°[A, S, D, G, V, U, Y] + F int [A, S, D, G, V, U, Y] . 

We use the following notation: V is the self consistent quark-gluon vertex, U is the self consistent 3-gluon vertex, and 
Y is the self consistent ghost-gluon vertex. These propagators and vertices are to be determined self-consistently from 
the equations of motion. We also define i> = i(T° + F nt ). We show $ graphically in Fig. HI We note that shifting 
the field introduces another 3-point vertex and that, in the figure, the intersection of three gluon lines at a small dot 
represents fi° fc + A t M? jkl . 



24 



e 



e-,e 



fVWWAA^ -| ^VWIAA^ 









FIG. 5: The 3-loop 3PI effective action. Wiggly lines represent gluons, solid lines are quarks and dotted lines are ghosts. 

The equations of motion are obtained from the stationarity of the action. There are 11 equations which are obtained 
by functionally differentiating with respect to the 11 functional arguments of the effective action: 



ST 

5X t 



0; Xie{il>,<$,A,ri,fj,S,D,G,U,V,Y}. 



(17) 



The equations obtained by varying with respect to {S, D, G, U, V, Y} can be solved simultaneously for the self con- 
sistent solutions which are functions of the field expectation values: S[i^, A, rj, 77], D[^, A, 77, 77], G[^,^,A,r},fj\, 
U[i/j, iji, A, 77, 77], V[il), ^, A, 77, 77], Y[^, ^, A, 77, 77]. We write this set of self consistent solutions: 



X = {S,D,G,U,V,Y}. 



(18) 



Substituting these self consistent solutions we obtain the resummed action, which depends only on the expectation 
values of the fields: 

f[7/;,7/U, 77 , 77] (19) 
= T[^, t/S, A, S[ip, 4), A, 77, 77], D[ip, A, 77, 77], G[i/j, A, 77, 77], Vty, A, 77, 77], U[i/>, A, 77, 77], Y[i/>, A, 77, 77]] . 

The equivalence of (p~6|) and (|T9|) at the exact level was shown in [27|. In the future we will write Y and T without 
their arguments. 



V. 'EXTERNAL' n-POINT FUNCTIONS 



We define some type (2) 'mixed' vertex functions using the same notation as (|T5|) : 



132 



sd£ 



e 



6 A3 ' ~ 132 6 A3 ' " 1J " <L4 3 
These vertices are shown in Fig. |6l 'External' gluon legs are distinguished by a star. 



A 



132 



(20) 



\ 




FIG. 6: Some 'mixed' vertices. 
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Some additional useful relations can be obtained from the identities: 

D^D^ = #12 > G 13 1 G32 = 8 12 j ^13 1 ^32 = #12 • 

Differentiating ([2Tj) with respect to A and using ([20|) gives: 

= ^ll'^l'32'^2'2 , 7"T = ^i]./ 3>i/ 32 / G 2 / 2 , f 7^ = ffll' A1/32/ #2'2 ■ 



The 'external' gluon propagator is defined as: 



(21) 



(22) 



z(^f )- 1 = 



-f[^,^,i4,?7,7/] . 



The 'external' self energy is extracted from the 'external' propagator using: 



(23) 



(24) 



We can derive an expression for the 'external' self energy as a function of the vertices in (|2Q|) by taking derivatives of 
the modified effective action and using the chain rule. We use the notation Xi to indicate one of the set of functional 
variables X := {S, D, G, V, U, Y} and Xi to indicate one of the set of self-consistent solutions X := {5, D, G, V, U, Y}. 
We obtain: 



5 2 T 



(25) 



5A 2 SAi 



x 



E 



ST 

5X~i 



S 2 X t 



<5 2 r 



SXj 
x5A 2 



£ 2 r 



SXiSXj 



x SAi SA 2 



The second term in this result is identically zero (see Eqn. (|T7|) ). The expression can be further simplified by using 
the set of equations obtained by differentiating the equations of motion: 



s 


\ 5T 




Ja 2 




X- 



= => 



s 2 r 



SXi5A 2 



x 



E 



5 2 T 



SXjSXi 



xSA 2 



Using this constraint (|25|) becomes: 



s 2 r 



SA 2 5Ai 



x 



E 



s 2 r 



SXrfAi 



SXj 
xSA 2 



(26) 



(27) 



Expanding the sum and using (22) we have: 

s 2 r 



S 2 T 



SA 2 SAi 



x 



SSs^SAi 
S 2 T 



(S , 33/A3/24 / ^4 / 4) 



S 2 T 



Using (p~6|) we obtain: 



345 



5 2 T 



X 5A 2 5V345^1 



6D346A1 

^ 

X 5A 2 Jy 3 45^i 



(D33/ 1^3/24' ^4'4) 



S 2 T 



5Y ? 



345 



x SA 2 



. <5 2 r 



<5S 3 4<5Ai 

<5 2 r 



(28) 



5 2 r 



^G34^Ai 



(G33/ 63/24' G4'4) 



^34^1 

£ 2 r 



<5 2 $ 

/ s 2 ^ 

_ • (6 , 33/A 3 /24 / ^'4 / 4) = A — 

x \ 00 / 413 

1 / £ 2 $ 
„ • (D33'fi 3 '24'£>4'4) = O + 2—— 

x z V oUM 



5G34M1 



J 2 $ 

x ( G ^e 3 , 24 ,G 4 , 4 ) = I Go - 



(SAS) 32 4 
413 

(G9G) 3 24 • 



413 



(29) 
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The last three terms in (|28|) do not contribute because the derivatives SU /SA, SV/SA, SY/SA correspond to effective 
4-point vertices, which are not part of our leading order calculation. We extract Ilf^ from (|24|) . (|28j) and (|29|) : 

nif = \m» 3A2 d 43 + \ (no + 2^) 4i3 • (z>nz>) 324 - (a - • (*as) 324 (so) 

-( eo -^4j 413 - (GeG)324 + ^- 
The terms in round brackets can be written collectively using the notation: 

* : = (* + • A ° : = ( A » - mx) ■ e ° ; = ( e ° - ■ < 311 

The result for the 'external' self energy in (|30|) is shown in Fig. [71 The open circles in the figure denote Qq, Aq and 
6q, and the solid dots are fl, A and 6. 



^12* — o 



5 2 <f> 
SA 2 SA 1 



FIG. 7: The 'external' self-energy. 
The last term in Fig. can be calculated using the expression for <l> shown in Fig. \S\ The result is shown in Fig. [8l 

FIG. 8: A contribution to the conductivity. 

In order to calculate the viscosity using the 3PI formalism we use Eqn. ([T4)h with the integrands for the pieces 
IIint[gluon]^| (K, Q), Hint [quark] s a s b (K, Q) and Hint [ghost] *| given by the 2nd, 3rd and 4th terms in (|30j) . which are 
shown in the 2nd, 3rd and 4th diagrams in Fig. The vertices ft, 6 and A, and the self consistent vertices U, V 
and Y, satisfy a set of coupled integral equations that resum the pinching and collinear singularities. In the next two 
sections we derive these integral equations. 

In most cases, the equations in this paper are easier to understand when represented diagramatically. From this 
point on, we will give most results only as diagrams. When equations are used, all indices are suppressed. 



VI. INTEGRAL EQUATIONS FOR 'MIXED' VERTEX FUNCTIONS 

In this section we derive the integral equations for the type (2) 'mixed' vertex functions A and that appear 
in the bubble diagrams in II ext (see Fig. [7)). As explained earlier, these equations are obtained by taking functional 
derivatives with respect to the field expectation values of the appropriate equations of motion ([17)) . 

The integral equation for the vertex Q is obtained from the equation: 



5 


'ST 




5A 


6D 


X- 



0. 



(32) 



The subscript X indicates that all self consistent solutions 
ward to show that this expression can be written: 



T8|) are substituted. Using ([T6|) and ([20]) it is straightfor 



SA 



' 5 2 <S> 




SD SXi 


x. 



(33) 
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where the summation indicates contributions from all of the terms in the set {S, D, G, V, U, Y}. The terms Xi G 
{U, V, Y} give no contribution, because the derivatives SXi/SA correspond to effective 4-point vertices, which are not 
part of the leading order calculation. Expanding the sum we have: 



n = n' n + 2 



SD 

Ja 



' <5 2 $ 






" S 2 <f> 




JG 


' S 2 <f> 




SDSD 


X- 


+ 2 TT 
SA 


SDSS 


X- 


5A 


SDSG 


X- 



(34) 



The integral equations for the vertices A and 9 are obtained in exactly the same way. The results are: 



A = An 



" S 2 § 




SS 


' S 2 <$> 




SG 


' s 2 § 




SGSD 


Jt- 


~ SA 


SGSS 


X. 


~ Ja 


SGSG 


X- 



(35) 



SD 


' s 2 $ 




ss 


- s 2 <$> 




SG 


' s 2 § 




Ja 


SSSD 


X- 


~ Ja 


ssss 


X- 


~ Ja 


SSSG 


X- 



To simplify the notation, we define the 4-point functions: 



M ss ■-- 



S 2 $ 
SSSS ' 



M SD ■-- 



l SDSS 



M 



DD 



5 2 $ 
SDSD 



(36) 



where the dots indicate that the definitions for the 4-point functions involving ghosts are defined like the ones with 
quark fields. From Fig. [5] it is clear that Msg = Mqs = 0. For clarity, we give one example in Fig. [9l with all indices 
written out explicitly. The legs on each side of the box will join to a pinching pair of propagators. 




FIG. 9: A 4-point vertex. 



Using this notation, equations (|34|) and (|35|) can be represented diagramatically as shown in Fig. \W\ The shaded 
boxes in the figure represent the 4-point vertices Mxy- We use the notation: Mdd = box with diagonal lines, Mds, 
Msd = box with hatched lines, Mdg, Mqd = light grey box, Mss = dark grey box, Mqg — hatched grey box. The 
open circles denote the vertices Q' , A f and O' . 




FIG. 10: Structure of the integral equations for the vertices Q, A and B. 



We can calculate the bare vertices directly from (|3Tj) . Since 5 2 &/5S SA = 5 2 &/5G5A = 0, we have Aq = Ao and 
©o = ©o- The vertex Q' is slightly more complicated. After many cancellations, the surviving terms are shown in 
Fig. HU We note that ft f contains Mq and U, but not ft. In Fig. [Til we have combined diagrams that correspond to 
permutations of external legs. The third and fifth diagrams on the right hand side have a loop insertion on the upper 
leg. These diagrams should each be drawn as two diagrams, with symmetry factor 1/2, one with the loop insertion 
on the upper leg and one with the loop insertion on the lower leg. 
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FIG. 11: The vertex tl' . 



Below we discuss in detail the calculation of Mdd- In Fig. [12] we show each term in and the corresponding 
contribution to Mjjjj. As in Fig. [TTJ we combine diagrams that correspond to permutations of external legs. The 
diagram on the right hand side of part (b) in Fig. [12] should be drawn as two diagrams, one with the dotted vertex at 
the top and one with the dotted vertex at the bottom. Similarly, the diagram in the second part of the right hand side 
of part (e) should be drawn as two diagrams with the triangular insertion at the top in one diagram and the bottom 
in the other. In the same way, the first and second diagrams in the right hand side of part (f), and the diagrams 
in right hand side of (g) and (i) should be drawn as two diagrams. The third diagram in the right hand side of (f) 
should be drawn as 4 diagrams. 




FIG. 12: Diagrammatic representation of Mdd. The crossed lines in parts (e), (h) and (j) pass over /under each other and do 
not intersect at a 4-point vertex. 
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In Fig. [13] we give the results for Msd and Mss- The corresponding results for Mqd and Mqg have exactly the 
same form as Msd and Mss respectively, with the quarks replaced by ghosts. 



FIG. 13: Diagrammatic representation of Msd and Mss- The crossed lines in the last two diagrams on the right hand side pass 
over/under each other, and do not intersect at a 4-point vertex. As in Fig. [12J we have combined diagrams that correspond to 
permutations of external legs 

The integral equations for the vertices Q, A and 6 are obtained by substituting the equations represented in Figs. 
[TT| n~2l and fT3l (and the corresponding equations for ghosts), into Eqns. (|34|) and (|35|) (shown in Fig. [TO]) . 



In this section we derive a set of integral equations for the self consistent vertex functions. We show in detail how 
each term is obtained for the integral equation for the vertex U. For V and Y we give only the final expression. 

The integral equation for the vertex U is obtained from the equation of motion ST/SU — 0. In Fig. [HI we list the 
terms in and the corresponding contributions to the integral equation. 




M ss 




VII. EQUATIONS OF MOTION FOR THE SELF CONSISTENT VERTICES 









3 



i 

>■■■■•' 



FIG. 14: Contributions to the integral equation for U. 



Combining and rearranging, we obtain the integral equation shown in Fig. [T5l As before, we have combined diagrams 
that correspond to permutations of external legs: the third diagram on the right hand side of Fig. [15] should be drawn 
as 3 diagrams, each with symmetry factor 1/2. 
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V-. 



FIG. 15: The integral equation for U. 

In exactly the same way, we can obtain integral equations for the vertices V and Y. The result for the vertex V is 
shown in Fig. [TB] The equation for the vertex Y has the same form with the quarks replaced by ghosts. 





FIG. 16: The integral equation for V. 

As a check of our results, we verify that the formalism does not double count the collinear divergences. We start 
by considering the diagrams on the right hand side of Fig. [T2j We substitute Uq into the diagram in part (b) using 
Fig. [TH The result is to change the sign of diagram (d), and to cancel the second diagram in part (e), the second 
and third diagrams in part (f), and the diagrams in part (g) and (i). This result is shown in Fig. [17] The diagrams 
that have been removed by the substitution are exactly cancelled because they are already contained in the diagrams 
shown in Fig. [T71 together with the integral equations for U, V and Y, as shown in Figs. [15] and [TBI 





WAAAAAA; *W# — » #A/ 

2S -X 



FIG. 17: A rearrangement of the result for Mdd. 
Similarly, using Fig. [TB]to remove Vb from Fig. [13] produces Fig. [18] 

M SD = 2 



M SS 




FIG. 18: A rearrangement of the results for Msd and Mss- 

Again, we see that the diagrams that have been removed are precisely those that are contained in Fig. [16] The 
cancellation works in exactly the same way for Mqd and Mgg- 



VIII. SCATTERING AND PRODUCTION PROCESSES 



The results of the previous sections can be summarized as follows. In order to calculate the leading order QCD 
shear viscosity using the 3PI formalism we use Eqn. ([T4|h with the different pieces of the integrand for the self energy 
given in (|3Q|) (and shown in Fig. [7]). The vertices ^, 6 and A and the self consistent vertices U, V and Y satisfy 
a set of coupled integral equations that resum the pinching and collinear singularities. These integral equations are 
produced naturally by the 3PI formalism and are shown in Figs. \W\ fTol and [TBI 
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We can demonstrate that the integral equations produced by the 3PI formalism are correct by showing that the 
kernels of the Q, and A equations have the form of the square of the sum of the amplitudes that correspond to all 
physical scattering and production processes. In this section we outline the strategy of the calculation. Some details 
are given in Appendix The basic steps are as follows. 

(1) We re-expand the equations shown in Fig. [lOj keeping all terms to 2-loop order. Equivalently, we keep all 
contributions to the 4-point functions, up to 1-loop order. The expansions of the shaded boxes are shown in Fig. [T3 
and Fig. [TBI We use Figs. [15] and [16] to expand the U, V and Y vertex functions, and propagators are expanded by 
inserting the 1-loop pieces of the self energy (see, for example, Figs. [20l and |2T|) . 

(2) We perform the summations over Keldysh indices. The calculation for the A vertex for QED is done in detail 
in [15|. The method used here is exactly the same. We can write the resulting set of equations in the form: 

V X (P)=V°(P) + Yl S xy JdKM xy (P,K)A r y et (K)V y (K)A a v dv (K); x £ {gl, gh, q} , (37) 

2/€{gl,gh,q} 

v# = n, v gh = e, ^ q = A, 

A g i = D , A gh = G, A q = S, 

Mgigi = Mdd , Mgigh = Mjjg 5 M g \ q = Mds •> etc, 

where 'gl' stands for gluon, 'gh' stands for ghost, and 'q' stands for quark. We have supressed all indices, except that 
we have explicitly written the momentum variables and the integral over the 4-momentum K. The subscript q refers 
to a quark of a distinct flavour, and the sum in (|37|) is over each flavour of quark. The 3-point vertex is retarded with 
respect to the middle leg. The 4-point function is: 

M xy (P, K) = M xy {13, P, K) + N B {K) (M^(5, P, K) - M xy (9, P, K) 

where the numerical arguments of the 4-point functions indicate Keldysh components. This notation is explained 
in Appendix A. The factor S xy is the symmetry factor of the diagram. For example: for x = gl we have <S g i g i = 
1/2; <S g igh = <S g iq = — 1, and Eqn. (|37|) becomes: 



Sl(P) = Sl°(P) + J dK (38) 

M DD (P,K)D ret (K)Q(K)D adv (K) - M DG (P, K)G ret (K)0(K)G adv (K) - M DS (P, K)S ret (K)A(K)S adv (K)] , 

which is the equation shown in the first line of Fig. [TO] 

(3) Since we are only interested in verifying that the correct matrix elements are produced, we use factors for bare 
propagators in all numerators, but use HTL self energies to regulate pinch singularities in denominators. We rewrite 
the pinching pairs of propagators: 

S L e (3 S a?/3> = ~# oc(3# oc'(3> ~ ~ ^ 7— , ^ret = -Tr £ re £ (K)) , (39) 

21mL re t(A) ^ 
p(K) 

rivet riadv 

^ Xr ~~ 9 9 2Imn^(K) ' ^ ^ ~ ~ 2lmUj et (K) ' ret ~ ^ ret ' 

-ip(K) = 1/(K 2 + iSign(k )e) - 1/(K 2 - iSign(k )e) . 

Note that we regulate the pinching singularities from the pair of ghost propagators and the pair of gluon propagators 
with the transverse part of the gluon polarization tensor. This is justified because of the fact that after all cancellations 
have been taken into account, only transverse gluons survive. For future use we define: 

u x : x e {gl, gh, q} -> n g i = n gh = u^ et , n q = t ret . 



(4) From ([6]) we need the real part of each V x and consequently, from (|37|) . we need to extract the real part of each 
4-point function M xy (P, K). Our method is related to the Cutkosky rules at finite temperature and is described in 
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[28(| . Terms with an even number of on shell or 'cut' propagators are real. There are no terms with zero cut lines. 
It is easy to show that terms with four cut lines do not contribute, because it is kinematically forbidden for three 
on-shell lines to meet at a vertex: 

p(±(P 1 ±P 2 ))p(P 1 )p(P 2 )=0. (40) 

The conclusion is that all terms must contain two cut lines. However, some terms containing two cut lines are 
identically zero. Since the leg momenta P and K are on-shell (because they connect to pinching pairs of propagators) 
we can use (|4Q|) to obtain: 

p{K-L)A™\L) - p(K - L)Prm x (L) , p(K - L) A a x dv (L) - p{K - L)Prm x {L) , 

(41) 

A r x et (P-K) Prm x (P-K), A a x dv (P-K) Prin^P - K) , 

where we have defined Prin^(P) = l/2(A£ e *(P) + A x dv (P)). The result is that all non-zero terms contain two cut 
lines that effectively divide the diagram into the product of two amplitudes. 

We note that for diagrams where all propagators carry different momenta, the procedure described above is 
perfectly straightforward. In diagrams where more than one propagator carries a given momentum, one must be 
careful to show that potentially dangerous terms that contain the square of a delta function do not appear. The 
disappearance of these unphysical terms is a well known result due to the KMS condition (29| . 

(5) For each of the vertices V x , we define a new vertex V x by contracting each pinching line by the corresponding 
external leg. To simplify the form of later results, we also divide by the HTL width of the line: 



n = &(P) = [ - g Xr ] ^L- , = S"(P) = , i> = A"(p) = Tr [f a] ^ . (42) 



We recall that the ghosts are unphysical degrees of freedom whose only role is to cancel the contributions from the 
unphysical gluon polarizations. From Fig. 0we see that it is the combination V^i — 2V g h that appears in the viscosity. 
As a consequence, we will look at the vertices: 

V q , %, V g = V g i-2V gh , (43) 



where the vertex Vq is obtained from V q by conjugation. To distinguish these vertices from those defined in (|42|) . 
we use the indices {a, 6, • • •} G {<?, g}, instead of {x,y, • • •} G {q, q, gl, gh}. Our goal is to rewrite the integral 
equations for the vertices V x (in Eqn. (|37|) ) in terms of the vertices V a , and to show that these equations have the 
correct form, with the kernels given by the square of the sum of the amplitudes that correspond to the relevant 2 — > 
2 scattering and production processes. 

(6) In order to obtain the traditional form of the matrix elements, we must label the momenta in a specific way. 
Each contribution to the 4-point functions M xy has the form of a cut 1-loop amplitude. Each amplitude depends on the 
two external momenta P and K, and one loop momentum variable that is integrated over. We can introduce a second 
momentum integration by adding a 4-dimensional delta function. We relabel these four momenta by the four variables 
{P, P2, Pi, P2} which are defined so that they correspond to the two external momenta, and the two momenta carried 
by the cut lines (as discussed in step (4) above, there are always two cut lines). In addition, we choose directions so that 
P+P2 = P1+P2, which means that the final expression will contain an overall factor j dLi J (IL2 S 4 (P+P2 — Pi — L2). 
In principle, there are 16 terms which correspond to the 2 4 possible choices for the signs of the 0-components of the 
momenta on the four on-shell lines. Since P is an external variable, we make the choice p° > 0, which leaves eight 
terms. Only three of these terms correspond to kinematically allowed 2^2 scattering and production processes. 
For each diagram we write one of these three terms by choosing Sign(p°) = Sign^) = Sign(Zj) = Sign^). The 
terms corresponding to the other two choices can be obtained by making the changes of variables: P2 <-> —Pi and 
P2 — P2. We define the notation: 

/ dP 2 j dLi I dL 2 5\P + P 2 -L 1 - P 2 ) p 25 ^1,^2) (44) 

J J J p erms 

= f dP 2 J dL x J dL 2 5 4 (P + P 2 -L 1 - La) (j(P, P 2 ;L U L 2 ) + f(P, -L i; -P 2 , L 2 ) + f(P, -L 2 ; L u -P 2 )j . 
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Using the notation defined in (|39|) . ([42]) and (|44|) we can rewrite (|37|) as: 
0(p o ) 2Imn a K = fl(po) 2Imn a T>° - 0(p o ) ^ / • [ [ 

verms L J l 2 (,\ 



(45) 



M a { ^ cd (P,P 2 ;L u L 2 ) • (^^(L^Laj^Pa))^^^) 



H(P 2 ) • (2tt) a 5\P + P 2 -L 1 - L 2 ) . 



The sum over (z) in (|45|) includes the pairs of amplitudes M • M produced from the cuts of all of the 4-point functions 
that appear in the expansion of the kernel of the integral equation. The indices {ra a , l bl z c , jd} take the values b or 
f, depending on whether the corresponding line is a boson or fermion. We define: 

N(ijlm) = ab(l, E P2 )em(i, E h )em(j, E h ) /em(m, E p ) , (46) 

em(6, E x ) = l + n b (E x ) , em(/, E x ) = 1 - rif (^) , ab(6, ^) = n b (E x ) , ab(/, ^) = n f (E x ) . 

In Fig. [19] we give two examples of the way in which cutting a 4-point function produces the product of two 
amplitudes. The dashed line indicates the cut propagators. The cut box graph gives the square of the t-channel. The 
cut triangle graph gives the product of the t-channel and the s-channel. 

P 2 P Lj L, p 





'Li L 2 - 










p 






p 




P2 




L 


,y L2 





P2 



P2 



P2 L 2 



L 2 P 2 



A, 2 P 2 X 



FIG. 19: Two examples of a cut 4-point funtion written as the product of two amplitudes. 



The goal is to show that the quantity in square brackets in (45) can be rewritten in the form: 

2 



V b (P 2 



V a Jl X Jl 2 f , r 

1 2 {bcd}e{gq 



M 



ab^-cd 



(P, P 2 \L\,L 2 \ 



N(mJ b i c j d ) 



V b (P 2 ). 



(47) 



The amplitude J\/[ ah ^ cd denotes a scattering amplitude for the process ab — > cd and the square is summed (not 
averaged) over the spins and colours of all states. As before, the subscripts q and q denote quarks and anti-quarks of 
distinct flavours, and the sum is over each flavour of quark. . The factor v a is equal to the number of spin x colour 
states for the external excitation, so that dividing by v a produces an average over initial states. The thermal factors 
give the correct combination of statistical emission and absorption factors. Since we have assumed Sign(po) > 0, we 
obtain the loss term. The choice Sign(po) < would produce the gain term. 



If we take a = g in (|45|) we obtain: 

Vg Jl ± Jl 2 

[M 99 " 99 N(bbbb) + M"^ qq ~N(bbf i)]V g (P 2 ) + M 9q ^ 9q N(bfbf)V g (P 2 ) + M 9q ^ 9q ^{bM)Vq{P 2 ) 
l&d A C\ (3 



(48) 
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Taking a = q in (jUJl we obtain: 



V b {P 2 ) = - 

M qq ^ qq N(f ff f)V q (P 2 ) + M g9 ^ g9 N(fbfb)V g (P 2 ) + [M qq ^ qq N(f f f f) + M qq ^"N(i fbb)]V r ? (P 2 ) 



(49) 
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The numerical subscripts in Eqns. (|48|) and (|49|) refer to quark flavours. Our results agree with the SU(3) results of 
[30j |. and with the results of [3]. Some details of the calculation of Eqn. (|48|) are given in Appendix O 

IX. CONCLUSIONS 

In this paper we have presented the first calculation of the complete leading order QCD shear viscosity using 
quantum field theory methods. We have demonstrated that the calculation can be organized naturally using the 
3PI effective action. The expression produced by the Kubo formula (Eqn. (fT3|) ) contains vertices which satisfy 
a set of coupled integral equations that resum the pinching and collinear singularities. These integral equations 
are produced naturally by the 3PI formalism, without the need for any kind of power counting arguments, and 
are shown in Figs. [TUl HU and HU We have verified that the integral equations produced by the 3PI formalism 
are correct, by showing that the kernels of the Q, 6 and A equations have the form of the square of the sum 
of the amplitudes that correspond to all physical scattering and production processes. In principle, the method 
developed in this paper should be generalizable to the calculation of transport coefficients at higher orders. 
Work in this direction is in progress. Our calculation provides a connection between nPI effective theories and 
kinetic theories, and supports the use of nPI effective theories as a method to study the equilibration of quantum fields. 



APPENDIX A: KELDYSH REPRESENTATION 



We use the closed time path formulation of real time statistical field theory [3l|, |32| which consists of a contour 
with two branches: one runs from minus infinity to infinity along the real axis, the other runs back from infinity to 
minus infinity just below the real axis (for reviews see, for example, [H, [H]). The closed time path contour results in 
a doubling of degrees of freedom. Physically, these extra contributions come from the additional processes that are 
present when the system interacts with a medium, instead of sitting in a vacuum. As a result of these extra degrees of 
freedom, n-point functions have a tensor stucture. Statistical field theory can be formulated in different bases, which 
produce different representations of these tensors. We will work in the Keldysh basis. In the discussion below, we 
use bi = 1 or 2 to denote indices in the 1-2 basis and ci to denote indices in the Keldysh basis, with a = 1 := r and 
Ci = 2 := a. The rotation from the 1-2 representation to the Keldysh representation is accomplished by using the 
transformation matrix: 



Keldysh-^- (1-2) 



1 

71 



(Al) 
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The vertices in the Keldysh representation are given by: 

In order to simplify the notation for the vertices, we replace each combination of the indices {r, a} by a single numerical 
index. In momentum space we write: 

r ClC2 - c "( Pl ,p 2 , • • - Pn ) = T(i, Pl , P2 , ■■■pn). (A3) 

We assign the choices of the variables c\c<i ■ • ■ c n to the variable i using the vector: 

where the symbol <g> indicates an outer product. For each n, the zth component of the vector corresponds to a list of 
variables that is assigned the number i. To simplify the notation we drop the subscripts and write a list like r^as 
as rra. For clarity, the results are listed below. 

3- point functions: rrr — > 1, arr — > 2, rar — > 3, aar — > 4, rra — > 5, ara — » 6, raa — » 7, aaa — > 8, 

4- point functions: rrrr — » 1, arrr — > 2, rarr — » 3, aarr — > 4, rrar — > 5, arar — > 6, raar — > 7, aaar — > 8, rrra — > 9, 
arra ^ rara H 5 aara — > 12, rraa — > 13, araa — > 14, raaa — > 15, aaaa — > 16. 

Summations over Keldysh indices can be done by hand, but the process is extremely tedious. In- 
stead, we use a Mathematica program. This program is described in detail in [34] and is available at 
www.brandonu.ca/physics/fugleberg/Research/Dick.html. The program can be used to calculate the integrand 
corresponding to any diagram (up to five external legs) in the Keldysh, RA or 1-2 basis. The user supplies input in 
the form of lists of momenta and vertices for each propagator and vertex. 



APPENDIX B: EXPANSION OF THE 'EXTERNAL' 2-POINT FUNCTION 

As a further check on our calculation, we can show explicitly that to 2-loop order the 'external' 2-point function in 
Eqn. (|24|) (shown in Fig. [7)) contains all terms that one would obtain from a straightforward Wick expansion, with 
the correct symmetry factors [35|. We start with the second term in Fig. [71 The bare vertex on the left hand side is 
f^Q as shown in Fig. [Til The dotted vertex on the right hand side is given by the first integral equation in Fig. \10\ 
and Figs. [T2land[T3l Term by term we obtain: 
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^%DDM DD DD% 
^%DDM DS SSA' 
^Q f DDM DG GGe 



FIG. 20: Some 2-loop contributions to II ext . 

As discussed before, we have simplified the figure by combining diagrams that correspond to permutations of external 
legs: the last diagram in the first line of Fig. [20l should be drawn as 2 diagrams, each with a symmetry factor of 1/2, 
with the loop insertion on the right and left hand side. The diagrams in the second and third lines are obtained by 
expanding self consistent propagators and inserting the first diagram in Fig. El and the 1-loop diagrams in Figs. [20l 
and EH Since we are working to 2-loop order, the full vertices U, V and Y in all 2-loop diagrams can be immediately 
replaced with the bare vertices, using Figs. [15] and [TBI 

Now we look at the third term in Fig. The dotted vertex on the right hand side is given by the second integral 
equation in Fig. [TUl and Fig. [T3l Term by term we obtain: 

\^oSSA f 



-AoSSMsdDD^ 
^A SSM ss SSA f 



FIG. 21: Some 2-loop contributions to II ext . 

The third diagram in Fig. Ogives the same diagrams as in Fig. [211 with the quarks replaced by ghosts. 

To get the full self energy to 2-loop order we combine all contributions. We need the set of diagrams in Figs. [20] 
and [21] (and the corresponding terms for ghosts). We also need the tadpole graph in Fig. [71 and the two loop graph 
obtained by expanding the self consistent propagator in the tadpole graph and re-inserting the tadpole. This produces 
the double scoop diagram shown in Fig. [22] Finally, we also need the sunset diagram in Fig. [8] The symmetry factors 
are all given explicitly on all diagrams. The resulting set of diagrams is the complete result for the 2-point function, 
at two loop order. 
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FIG. 22: The double scoop contribution to II' 



APPENDIX C: THE V g EQUATION 

In this Appendix we give some details of the calculation of Eqn. (|48j) . We follow the steps outlined in section [Villi 

(1) We iterate the integral equations and keep all terms up to 2-loop order, or contributions to the 4-point functions 
up to 1-loop order. These graphs can be divided into types based on their topologies. We call them box, bubble, 
triangle, cross, loly, tent and fish graphs. They are shown in Figs. [23j EH [25j [26] and [271 respectively. Each graph 
carries a numerical factor that is not included in the figure. Graphs that are labeled with the same number, such as 
(5) and (5') in Fig. [23j give the same result but are drawn separately so that it is easier to see that all contributions 
are included. The numerical factors are listed below, in the same order as the graphs in the corresponding figure. For 
example: the first line in (|C1|) means that the numerical factor for the first box graph is 1, the factor for the second 
box graph is -2, etc. 



factor[box] = {1, -2, -2, -2, -2, -2, -2, -2, -2, 4, 4, -2} , 
factor[bub] = j^, -1, -1, -1, -2, -2, -2, 2, 2 j , 
factor [tri] = {2, -4, -4, -4, -2, -2, -8, -8, -4, -4} , 



(CI) 



factor [cross] 



factor[loly, tent, fish] = |l,2, 

In Figs. [23l EH [25] and [271 the dashed lines indicate the two internal propagators that are cut, as discussed in section 
IVIIIi For the crossed-box graphs, there are two possible combinations of on shell internal lines which correspond to a 
horizontal and a vertical cut. In order to simplify Fig. [26j we do not draw the dashed lines that correspond to these 
cuts. 



(i) 



(5) 




(2) 



(6) 



(3) 



(7) 



«WV» 'WV. 
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FIG. 23: Box graphs that contribute to the integral equation for V g . 
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(1) 




(2) 



(2') 



(3) 




(4) 



(5) (5') — (6) 



(7) 



FIG. 24: Bubble graphs that contribute to the integral equation for V g . 



(i) 
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FIG. 25: Triangle graphs that contribute to the integral equation for V g . 



(i) 




(2) 




(3) 



(4) 




(5) 




(6) 




(7) 




FIG. 26: Crossed-box graphs that contribute to the integral equation for V g . 



FIG. 27: Loly, tent and fish graphs that contribute to the integral equation for V g . 

(2) We calculate the Keldysh structure for each diagram. In order to do this, we look at the corresponding diagram 
where all lines are scalars, in the sense that they have no Dirac or Lorentz structure, but carry the appropriate boson or 
fermion thermal distribution functions. The summations over Keldysh indices are done using a Mathematica program. 
The program is described in [34] and is available at www.brandonu.ca/physics/fugleberg/Research/Dick.html. The 
program can be used to calculate the integrand corresponding to any diagram (up to five external legs) in the Keldysh, 
RA or 1-2 basis. The user supplies input in the form of lists of momenta and vertices for each propagator and vertex. 
Several examples of this part of the calculation are worked out in detail in [15|]. We extract the overall phase space 
factor: 

T= lljl 2 := {2*Y2E h (2^2^' (C2) 

and define the mandelstam variables: 



s=(L 1 +L 2 ) 2 ; t = {P-L 1 f- u=(P-L 2 ) 2 . (C3) 



21 



The results are listed below. The notation 'den' is a reminder that the contributions from the numerators, produced 
by correctly including the appropriate Dirac and Lorentz structure, are not yet included. 

den[box(l)] = ^N(bbbb) • 1 ; z G {1, 5, 6, 7, 9, 10} , (C4) 
den[box(2)] = ^N(bbfF) • ^ , den[box(3)] = den[box(4)] = den[box(8)] = - ^N(bfbf) • ^ , 

den[bub(l)l = J^Nfbbbb) • \ ; z G {1, 2, 5, 7} , 

s z 

den[bub(3)l = den[bub(6)l = ^N(bbff) • \ , den[bub(4)l = - ^N(bfbf) • \ , 
den[tri(l)] = J^N(bbbb) • ; i E {1, 3, 6, 7, 8, 9} , 

den[tri(2)] = J r N(bbff) • - , den[tri(4)] = - J r N(bfbf) • - , den[tri(5)] = - J r N(bfbf) • - , 

st st st 

den[cross(l)] = J^N(bbbb) • ^- ; i e {1,3,5, 6, 7} , 

den[cross(2)] = ^N(bbff) • — , den[cross(4)] = - ^N(bfbf) • — , 

tu tu 

den[loly] = Ti N(bbbb) • ^ , den [tent] = Ti N(bbbb) • j , den[fish] = -J^N(bbbb) . 

The results for the numerators are given in (|C5|) and are listed in the same order as the diagrams in the corresponding 
figure. For example, the first line in (|C5|) gives the numerators for the box graphs labeled (1), (2), etc. The crossed-box 
graphs produce two contributions each, because there is a horizontal and a vertical cut. 

num[box] = { 8 - 25suj C 2 A , -StuC 2 F d F , -SstC 2 F d F , 4 (4s 2 + 3us + 4u 2 ) C A C F d F , (C5) 

- 2 (t 2 + lOsu) C\, 2t 2 C 2 A , 2t 2 C\, -16suC 3 A C F d F , 4 ^ S — + ^ ^ C 2 A , 2t 2 C\ } , 
/69s 2 \ 

num[bub] = { 8 ( — 2htu\C\, 2 (-s 2 - lOtu) C\, 4 (At 2 + 3ut + 4u 2 ) C A C F d F , -SstC 2 F d F ,2s 2 C 2 A , 

ft 2 u 2 \ 
-2tuC A C F d F , 4f- + — )C 2 A } , 

num[tri] = ^8 (lbsu- C\, 4u 2 C A C F d F , u 2 C 2 A ,0, 4s 2 C A C F d F , 0, -s 2 C 2 A , t 2 C 2 A , -u 2 C\^ , 

num[cross - horz] = |4 (^30tu - ^-^j C 2 A , 0, 0, 4*x 2 C A C F d F , -u 2 C 2 A , u 2 C 2 A , oj , 

„u m[ c r „ ss - ve rt| = {4 [mu - 2£) C5.0.0, 4rtWW F . -^cj, AS. A3} • 



num[loly /tent /fish] = {324, z s, 324 z 5, -864 C^} 



The last step is to sum all contributions and show that: 

factor[z]num[z]den[z] = T(M 99 ^ 99 N(bbbb) + M"^ qq ~ N(bbff) + M gq ^ gq N(bfbf ) + M 9q ^ 9q N(bfbf )) . (C6) 

(i) 

The sum over (z) is over all of the graphs in Figs. [23l EH [25J, [26] and [27] We use the fact that the definitions of the 
internal momenta L\ and L2 can always be reversed, or equivalently, that each term can be written in a symmetric 
form by interchanging t and u. Note that when we insert the result for |^ a6 ^ cd | into Eqn. (|45|) . we must introduce 
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an extra factor 1/2 if the final states c and d are not the same, to avoid double counting this contribution. These 
factors of two are indicated in square brackets in (|C7|) . The results are: 

|M—f = 16^i(3-|-£J-g), (C7) 
= [I]l6d, C F (<7, (1 + |) - g 4)) , 

i-^^r = i^™t = [^] le^c, (^ c-i, + f )) , 

in agreement with (|48|) . 
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